Vortex Nucleation in a Stirred Bose-Einstein Condensate 
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We studied the nucleation of vortices in a Bose-Einstein condensate stirred by a laser beam. We 
observed the vortex cores using time-of-flight absorption imaging. By varying the size of the stirrer, 
we observed either discrete resonances or a broad response as a function of the frequency of the 
stirrer's motion. Stirring beams small compared to the condensate size generated vortices below 
the critical rotation frequency for the nucleation of surface modes, suggesting a local mechanism of 
generation. In addition, we observed the centrifugal distortion of the condensate due to the rotating 
vortex lattice and found evidence for bent vortices. 

PACS 03.75.Fi, 67.40.Vs, 32.80.Pj 



Dissipation and turbulence in superfluid flow often in- 
volves the creation and subsequent motion of quantized 
vortices Since vortices are topological defects they 
may only be created in pairs, or can enter a system indi- 
vidually from its boundary. The nucleation process has 
been a subject of much theoretical interest 1^. Exper- 
iments with Bose-Einstein condensates (BEC) confined 
in atom traps are well suited to test theories of nucle- 
ation because the boundary of the condensate is well 
controlled, and vortices can be directly imaged. 

Vortex nucleation usually involves dynamical instabil- 
ities and superfluid turbulence Exceptions are the 
direct coupling between ground and vortex states for 
small condensates |^ and the engineering of states of 
quantized circulation by manipulating the phase of the 
wavefunction ^-|^. Turbulent flow can be created by 
perturbing the system with a time-dependent boundary, 
for example by a small laser beam [Qjsj, or by a rotating 
trap anisotropy |9-11]. The resulting vortices have been 
directly imaged l^-lCj. 

In a rotating frame, vortices are energetically favored 
above a critical rotational frequency Q.c Pp^]. Current 
theories suggest that the nucleation of vortices requires 
the creation of surface waves which "collapse" into a vor- 
tex state p^-p7|. This corresponds to the disappear- 
ance of the energy barrier for a vortex to enter the cloud 
|p|, P"4|Jl^ . For surface excitations with angular momen- 
tum ITi and frequency oji, the minimum rotational fre- 
quency is 

17s = min/(a;//0. (1) 

This corresponds to a Landau critical velocity for sur- 
face waves Vc = ^sRtFi where i?TF is the Thomas-Fermi 
radius I psf . 

Our main result is that vortices can be created below 
^Is by using a small, localized stirring beam. This in- 
dicates that the current surface mode analysis of vortex 
nucleation is incomplete. Recent experiments [^,^ ob- 
served vortices only at much higher frequencies, giving 
rise to a variety of theoretical models [0,^ . One expla- 
nation is that those experiments excited only the I = 2 



mode, which requires a higher drive frequency than Eqn. 
(|l|) [ p^ , pl|j2^ ]. We verified this prediction by stirring 
the condensate with anisotropics of different symmetries 
(Z — 2,3,4) and observing distinct resonance frequencies 
for vortex formation. If the condensate is stirred with a 
small beam, one would expect to couple to many modes I. 
Surprisingly, the small stirrer did not excite resonances, 
but could nonetheless generate vortices as effectively as 
a resonant drive. 

Our method of vortex generation has been outlined in 
previous work We start with nearly pure BECs (> 
90% condensate fraction) of up to 5x10^ sodium atoms 
in a cylindrically shaped loffe-Pritchard magnetic trap 
with radial and axial frequencies of cj^ = 27rx 83 Hz and 
LUz = 27rx 20 Hz, respectively. A radio frequency "shield" 
limited the magnetic trap depth to 50 kHz (2.3 /iK). The 
condensate chemical potential and peak density were 300 
nK and 4 x 10^'' cm^"^, respectively, corresponding to a 
healing length ^ ~ 0.2 /im. 

Vortices were generated by rotating the condensate 
along its long axis with a scanning blue-detuned laser 
beam (532 nm) |23|. We explored different beam waists 
between 5 and 25 /im. For the tightest focus, the peak 
optical dipole potential was 620 nK. We used scan radii 
as large as the Thomas-Fermi radius -Rtf, which var- 
ied from 27-30 /im. Multiple beam patterns (formed by 
rapidly scanning the laser beam from 1.5-10 kHz) were 
used to create vortices. The laser beam was left on dur- 
ing evaporation to damp out dipole motion of the con- 
densate. Immediately after producing a condensate we 
began the rotation for times of up to 500 ms, generating a 
vortex tangle. The laser beams were then instantly shut 
off and the cloud equilibrated for 500 ms, during which 
time the vortices crystallized into an Abrikosov lattice as 
shown in Fig. ^ and detailed in previous work [0. For 
small numbers of vortices the gas did not fully settle into 
a regular lattice before imaging. 

The vortex cores were observed using resonant absorp- 
tion imaging after 41 ms of ballistic expansion, which 
magnified them by 20 from their size ^ in the trap. As in 
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our previous work we imaged a 50-150 /im slice of atoms 
in the center of the cloud using spatially selective optical 
pumping into the F = 2 to F = 3 cycling transition pOl . 
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FIG. 1. Discrete resonances in vortex nucleation. The 
number of vortices created by multi-point patterns is shown. 
The condensate radius was Rtp = 28 ^im. Each data point 
is the average of three measurements. The arrows below 
the graph show the positions of the surface mode resonances 
(jji = u)r/Vl- The stirring times were 100 ms for the 2- and 

3- point data, and 300 ms for the 4-point data. Inset shows 
2-, 3-, and 4-point dipole potentials produced by a 25 fim waist 
laser beam imaged onto the CCD camera. The separation of 
the beams from the center is 25 pm for the 2-point pattern 
and 55 ^m for the 3- and 4-point patterns. The laser power 
per spot was 0.35 mW, 0.18 mW, and 0.15 mW for the 2-,3- 
and 4-point data, respectively. 

By varying the stirring parameters we explored differ- 
ent mechanisms for vortex nucleation. A large stirrer, 
with a beam waist comparable to the Thomas-Fermi ra- 
dius showed enhanced vortex generation at discrete fre- 
quencies. Fig. 1^ shows the number of vortices versus 
the frequency of rotation of the laser beam using 2-, 3- 
and 4-point patterns. The total laser beam power cor- 
responded to an optical dipole potential between 60 nK 
and 240 nK. The resonances were close to the frequencies 
of excitation of I =2, 3 and 4 surface modes (w; = ujr/ Vl) 
A second, higher resonance appeared in the 3- and 

4- point data, presumably due to additional coupling to 
the quadrupole (/ = 2) mode. The shift in the position 
of the peaks from the frequencies uji of elementary exci- 
tations is probably due to the presence of the stirrer and 
the vortices. 

Our results clearly show discrete resonances in the nu- 
cleation rate of vortices which depend on the geometry of 
the rotating perturbation. This confirms the role of dis- 
crete surface modes in vortex formation. A dependence 
on the symmetry of the stirrer (1-point versus 2-point) 
has also been explored in Paris [g2|. For longer stirring 
times and higher laser powers the condensate accommo- 
dated more vortices at all frequencies, and the resonances 
became less pronounced. 

A stirrer much smaller than the condensate size could 
generate vortices very rapidly ~ more than 100 vortices 



were created in 100 ms of rotation. Fig. ga shows the 
number of vortices produced using a 2-point pattern with 
a scan radius close to Rtf for various stirring times. 
Above 300 ms the angular momentum of the cloud ap- 
peared to saturate and even decreased, accompanied by 
visible heating of the cloud. 
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FIG. 2. Nonresonant nucleation using a small stirrer, (a) 
Average number of vortices created for different stirring times 
using a 2-point pattern positioned at the edge of the conden- 
sate. The beam waist, total power, and separation are 5.3 fim, 
0.16 mW, and 54 fim, respectively, (b) Effective lattice ro- 
tation frequency as a function of stir frequency. The lines 
in both graphs indicate the predictions of different models 
described in the text. 

The maximum number of vortices was roughly propor- 
tional to the stirring frequency. For a lattice rotating at 
frequency O, the quantized vortex lines are distributed 
with a uniform area density Uv = —, where k = h/M is 
the quantum of circulation and M the atomic mass [ p^ . 
Therefore, the number of vortices at a given rotation fre- 
quency should be 



N„ 



27rR^n/K 



(2) 



in a condensate of radius R. The straight line in Fig. ga 
assumes R — Rtf and that the lattice has equilibrated 
with the drive. In contrast to the large stirrer no res- 
onances were visible even when the number of vortices 
had not yet saturated. This suggests a different mecha- 
nism for vortex nucleation. Further evidence for this was 
obtained from the frequency and spatial dependences. 

For our experimental conditions, Eqn. |^) yields — 
0.25w,. = 21 Hz [|4). With the small stirrer, we ob- 
served vortices at frequencies as low as 7 Hz. Below 
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this frequency the velocity of the stirrer was not much 
larger than the residual dipole oscillation of the conden- 
sate. The rotational frequency below which a rectilinear 
vortex in the condensate center is energetically favored 
is 7 Hz ||_ 

In FigTH we varied the radius of the 2-point scan. The 
frequency of rotation was chosen to keep the linear ve- 
locity of the laser beam constant. Vortices could be gen- 
erated over a broad range of radii. The maximum num- 
ber was obtained at intermediate radii rather than the 
Thomas-Fermi surface, providing further evidence that 
surface excitations of the unperturbed condensate are not 
the dominant nucleation mechanism. The observed ra- 
dial dependence makes it unlikely that the thermal cloud 
plays a crucial role in the vortex nucleation, since at very 
low temperatures its maximum density occurs at the sur- 
face of the condensate. Indeed, we observed that finite 
temperatures increased the damping and led to fewer vor- 
tices. 
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FIG. 3. Spatial dependence of vortex nucleation. (a) 
Shown is the number of vortices as a function of the stir- 
ring radius for a constant stirring velocity (0.2 c). At the 
Thomas-Fermi radius, the angular frequency of the stirrer was 
17.5 Hz. The error bars reflect the scatter of 3 measurements 
at the same radius, (b) Axial phase contrast images showing 
the pure condensate with the aligned stirrer (1), the conden- 
sate during (2-7) and after (8-10) stirring. 

The nucleation mechanism for small stirrers may be 
related to our earlier experiments on the onset of dissi- 
pation in stirred condensates, where we observed a drag 
force at velocities above ~ 0.1c, where c is the speed of 
sound at the condensate center The friction with 

the moving stirrer causes an asymmetry in the density 
profile in front of and behind the laser beam. This has 
been directly imaged for linear motion . Similar flow 
field effects can be observed in Fig. |^, where they are 
clearly linked to the formation of vortices. Vortex pairs 



are predicted to arise from linear stirring |2^,^. When 
the laser beam moves in a circle, co-rotating vortices will 
be favored, whereas counter-rotating vortices will be ex- 
pelled from the system. 

For an object smaller than the healing length ^, the 
critical velocity for vortex formation occurs at the Lan- 
dau value given by Eqn. (|l|). For objects larger than ^ 
such as our laser beam, hydrodynamic flow around the 
stirrer can reduce the critical velocity relative to Eqn. 
(|l|) [ p6[ , which may explain our observation of vortices 
below 21 Hz. 
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FIG. 4. Centrifugal effects on a rotating condensate. 
Shown is the observed radial size of the condensate versus 
vortex number. The condensate size was determined from a 
parabolic fit to the time-of-flight distribution after 41.5 ms. 
(b,c) Comparison of a non-rotating condensate (i?tof ~ 570 
fj,m) and one with 160 vortices. Radial in-situ phase contrast 
imaging of a condensate at rest (d) and in rotation at a fre- 
quency of 75 Hz (e), showing the modified aspect ratio. The 
length of the condensate in (d) was 200 nm. 

In our earlier work we observed that the size of a 
condensate with vortices exceeded the size of the non- 
rotating condensate [|l0| . Here we study these centrifugal 
distortions quantitatively (also note (ll)). Fig. ^ shows 
the enhancement of the cloud size R in time-of-ffight by 
up to 25% due to additional rotational kinetic energy. 
The rotating condensate experiences a centrifugal po- 
tential — (l/2)M(r2r)^, which leads to an effective radial 
trapping frequency of ui^ — 17^. For a constant number 
of atoms TV, this increases the Thomas Fermi radius 



Rtf 



Rq 



(3) 



(1_ (11)2)3/10 

and reduces the mean-field interaction energy Emt- The 
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total release energy of the gas is then E = Eint + ^IcS^'^ + 
E^oNy. The 2nd term accounts for the rigid rotation of 
the lattice while the 3rd term is a quantum correction 
due to the kinetic energy of the cores, and is negligi- 
ble for large N^. The effective moment of inertia of the 
condensate is Ics = 2/7AIR^pN. The energy per unit 
length of a single vortex is Eyo — nnh'^ /M In [0.27 d/S,], 

— 1/2 

where d — Uy is the mean distance between vortices 
p8[ . When we average Eyo over the Thomas-Fermi dis- 
tribution we predict a 30% increase in E for 120 vortices, 
whereas the observed increase in is about 50%. This 
discrepancy is probably due to our selection of the central 
slice of the cigar, where the rotational energy is the high- 
est, and due to the unobserved axial expansion, which 
may depend on the angular momentum of the cloud. 

If we account for centrifugal effects and combine Eqns. 
and (|^) , we expect a divergence of the number of vor- 
tices near the trap frequency (dashed line in Fig. |2^). 
The deviation of the data from this line suggests that 
the condensate did not fully equilibrate with the rotat- 
ing drive. Taking into account the critical velocity for 
vortex nucleation, Vc — 0.1c psf , we expect the maxi- 
mum rotation frequency of the lattice to be — Vc/Rtf 
where is the frequency of the stirrer moving at ra- 
dius i?TF- Using the measured number of vortices, we 
can invert Eqns. (||) and (||) to derive the lattice rota- 
tion frequency, which is shown in Fig. along with 
the expected value fig — 27r x 6 Hz assuming a constant 
A''. The discrepancy can be partly attributed to loss of 
atom number due to heating by the stirrer, which was 
up to 30%. We can also derive from these equations the 
flow velocity at the edge of the condensate. For a lat- 
tice with 144 vortices, this velocity exceeded the speed 
of sound at the condensate center by 40%, in contrast to 
a recent suggestion that supersonic rotation speeds are 
unattainable p3| . 

At low rotational velocities, vortices should not be rec- 
tilinear as assumed in many theoretical calculations but 
bent 30|. Such bent vortices should have lower visibil- 
ity in our images due to the line of sight integration across 
the optically pumped condensate slice. Fig. || shows sev- 
eral examples. Some appear as vortex lattices with tilted 
vortex cores. Other images show structures reminiscent 
of half rings and coiled vortices. However, since these 
are time-of-flight images it is not obvious how some of 
the observed features are related to spatial structures in 
the trapped condensate. 

In conclusion, we have studied the nucleation of vor- 
tices in condensates stirred by rotating laser beams and 
identified two different nucleation mechanisms, surface 
modes for large stirrers, as well as local turbulence for 
small stirrers. So far, we have mainly concentrated on 
pure condensates, but it will be intriguing to study the 
vortex phase diagram at finite temperatures and 
the role of the thermal cloud in the decay of vortices pSl . 
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FIG. 5. Three-dimensional structure of vortices. Shown 
are several examples of time-of-flight pictures of condensates 
at low rotational frequencies, where "smeared-out" vortex 
cores and elongated features were observed. The radius of 
the condensate was 510 fim. 
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